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We study the time evolution of a superposition of product states of two dressed atoms in a 
spherical cavity in the situations of an arbitrarily large cavity (free space) and of a small one. In 
the large-cavity case, the system dissipates, whereas, for the small cavity, the system evolves in an 
oscillating way and never completely decays. We verify that the von Neumann entropy for such a 
system does not depend on time, nor on the size of the cavity. 



I. INTRODUCTION 

Stability is a main characteristic of quantum mechan- 
ical systems in absence of interaction. When interaction 
with an environment is introduced to such systems, they 
tend to dissipate. A material body, for instance, an ex- 
cited atom or molecule, or an excited nucleon, changes of 
state in reason of its interaction with the environment. 
The nature of the destabilization mechanism is in general 
model dependent and approximate. An account on the 
subject, in particular applied to the study of the Brow- 
nian motion, can be found for instance in Refs. [l|, 0. 
However, stability (or not) of quantum mechanical sys- 
tems is not due only to the absence (or presence) of in- 
teraction. For example, the behavior of atoms confined 
in small cavities is completely different from the behav- 
ior of an atom in free space or in a large cavity. In the 
first case, the decay process is inhibited by the presence 
of boundaries, a fact that has been pointed out since a 
long time ago in the literature , while in the second 
case it completely decays after a sufficiently long elapsed 
time. 

This phenomenon of inhibition of decay and related 
aspects have also been investigated in [mlOl l using a 
"dressed state" formalism introduced in [ll|. With this 
formalism one recovers the experimental observation that 
excited states of atoms in sufficiently small cavities are 
stable. In [^,0], formulas are obtained for the probability 
of an atom to remain excited for an infinitely long time, 
provided it is placed in a cavity of appropriate size. For 
an emission frequency in the visible red, the size of such 
cavity is in g ood agreement with experimental observa- 
tions [13, Il3| . The dressed state formalism accounts for 
the fact that, for instance, a charged physical particle 
is always coupled to the gauge field; in other words, it 
is always "dressed" by a cloud of field quanta. In gen- 
eral, for a system of matter particles, the idea is that the 
particles are coupled to an environment, which is usu- 
ally modeled in two equivalent ways: either to represent 
it by a free field, as was done in Refs. [J 0|: or to con- 
sider the environment as a reservoir composed of a large 
number of noninteracting harmonic oscillators (see, for 
instance, (l^ - [T7| ). In both cases, exactly the same type 
of argument given above in the case of a charged particle 
applies to such systems. We may speak of the "dressing" 
of the set of particles by the ensemble of the harmonic 
modes of the environment. It should be noticed that 



our dressed states are not the same as those employed in 
optics and in the realm of general physics usually associ- 
ated to normal coordinates [3 [l3 ■ Our dressed states 
are given in terms of our dressed coordinates and can be 
viewed as a rigorous version of these dressing procedures, 
in the context of the model employed here (see Eqs. 
and pop in the next section). 

In the present paper we study the time evolution of a 
two-atom dressed state. This generalizes a previous work 
dealing with the simpler situation of a superposition of 
states of just one atom Our approach to this prob- 
lem makes use of the above mentioned concept of dressed 
states. We will consider our system as consisting of two 
atoms, each one of them interacting independently with 
an environment provided by the harmonic modes of a 
field. The whole system is supposed to reside in a spher- 
ical cavity of radius R. We take it as a bipartite system, 
each subsystem consisting of one of the dressed atoms. 
We will consider a superposition of two kinds of states: 
either all entities (both atoms and the field modes) are 
in their ground states, or just one of the atoms lies in its 
first excited state, the other one and all the field modes 
being in their ground states. The analysis of the density 
matrix of the system leads to the time evolution of the su- 
perposed states. The computation of the von Neumann 
entropy leads to the result that it remains unchanged 
as the system evolves, for a cavity of any size. We find 
rather contrasting behaviors for the time evolution of the 
system for a very large cavity (free space, i? — >■ oo) or for 
a small cavity. In the first case, as time goes on, the 
system dissipates completely, while for a small cavity the 
departure from the idempotency of the density matrix 
exhibits an oscillatory behavior, never reaching zero. 

The dressing formalism for just one atom inside a cav- 
ity is briefly reviewed in Section 2 in order to establish 
basic notation and formulas for the time evolution of the 
states. In Section 3, the formalism is generalized for the 
two-atom system and describe the evolution of its den- 
sity matrix, either in the case of a very large cavity (with 
infinite radius, that is, free space) or of a small cavity. In 
Section 4, we present our conclusions. 



II. DRESSING A SINGLE ATOM 

Let us briefly recall here some results from the analy- 
sis of previous works for the simpler situation of just one 
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atom, dressed by its interaction with the environment 
field. We shall thus consider an atom in the harmonic 
approximation, linearly coupled to an environment mod- 
eled by the infinite set of harmonic modes of a scalar 
field, inside a spherical cavity. A nonperturbative study 
of the time evolution of such a system is implemented 
by means of dressed states and dressed coordinates. We 
present, in this section, a short review of this formalism, 
for details see [ll| or [g^]. We consider an atom labeled 
A, having bare frequency uj\, linearly coupled to a field 
described by N (— )■ cx>) oscillators, with frequencies oj^, 
k — 1,2,...,N. The whole system is contained in a 
perfectly refiecting spherical cavity of radius R, the free 
space corresponding to the limit i? — > oo. Denoting by 
q\{t) {p\{t)) and qk{t) {pk{t)) the coordinates (momenta) 
associated with the atom and the field oscillators, respec- 
tively, the Hamiltonian of the system is taken as 
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where rjx is a constant and the limit N ^ oo will be 
understood later on. The Hamiltonian ([1]) can be turned 
to principal axis by means of a point transformation 
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performed by an orthonormal matrix T ~ (t^^^), where 
fi = (A, {k}) , k = 1,2,...,N, and rx = 0,...,N. The 
subscripts = A and ^ = k refer respectively to the 
atom and the harmonic modes of the field and rx refers 
to the normal modes. In terms of normal momenta and 
coordinates, the transformed Hamiltonian reads 
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where the U,r^ 's are the normal frequencies corresponding 
to the collective stable oscillation modes of the coupled 
system. 

Using the coordinate transformation, Ec^. in 
the equations of motion derived from the Hamiltonian 
Eq. ([T]), and explicitly making use of the normalization 

condition Y^^=a (^Ji^)^ = 1, we get 
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The right-hand side of equation ([5]) diverges in the limit 
iV — > oo. Defining the counterterm Suj'^ = Nrj\, it can be 



rewritten in the form 
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Eq. ([6]) has N + 1 solutions, correspondin g to the N +1 
normal collective modes. It can be shown [131 that if 
> Suj'^, all possible solutions for O^ are positive, phys- 
ically meaning that the system oscillates harmonically in 
all its modes. On the other hand, when < Suj'^, one 
of the solutions is negative and so no stationary configu- 
ration is allowed. 

Therefore, we just consider the situation in which all 
normal modes are harmonic, which corresponds to the 
first case above, > Suj'^, and define the renormalized 
frequency 



= j™^ (^A - ^^a) , 

following the pioneering work of Ref. j21j 
N ^ oo, equation © becomes 
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In the limit 
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We see that, in this limit, the above procedure is exactly 
the analogous of mass renormalization in quantum field 
theory: the addition of a counterterm —Nrjj^qJ^ {N 00) 
allows one to compensate the infinity of oj^ in such a way 
as to leave a finite, physically meaningful, renormalized 
frequency ujx- 

To proceed, we take the constant rjx as 
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where Aw is the interval between two neighboring field 
frequencies and g is the coupling constant with dimension 
of frequency. The environment frequencies uj^ can be 
written in the form 



TTC 



= fc — , fc = 1, 2, . . . , 

and, so. Aw — ttc/ R. Then, using the identity 
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Eq. ([SJ can be written in closed form 

'Rn\ _ n c 
~ J ^ 2^^ m 

The elements of the transformation matrix, turning the 
atom-field system to principal axis, are obtained in terms 
of the physically meaningful quantities flj-x ^^'^ wa after 
some long but straighforward manipulations [TT| , 
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The eigenstates of the system atoni(A)~field, 
\l\ThihT--)j a-re represented by the normahzed 
eigenfunctions in terms of the normal coordinates 
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xr^e~*^=A('=^+^)''=^*, (14) 
where Hi^^ stands for the ^^^-th Hermite polynomial and 

is the normalized vacuum eigenfunction, J\f\ being the 
normalization factor. 

Next, dressed coordinates q'^ and {<z^.} for the dressed 
atom and the dressed field, respectively, are introduced, 
defined by 



(15) 



{lox, LUk}. In terms of the dressed coordi- 



where u) 

nates, we define for a fixed instant, t — 0, dressed states, 
Ki, K2t ■ ■ ■) by means of the complete orthonormal set 
of functions UM 
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where, as before, /i labels collectively the dressed atom A 
and the field modes k — 1,2,3,..., that is, q'^ = q'^, {I'k}- 
The ground state Fq in the above equation is the same as 
in Eq. HH). The invariance of the ground state is due to 
our definition of dressed coordinates given by Eq. (1151) . 
Notice that the introduction of the dressed coordinates 
implies, differently from the bare vacuum, the stability 
of the dressed vacuum state since, by construction, it is 
identical to the ground state of the interacting Hamil- 
tonian in terms of normal coordinates. Each function 
V'kaki- - (?') describes a state in which the dressed oscilla- 
tor q'^ is in its K^-th excited state. 

The particular dressed state |r5'(0)) at t — 0, repre- 
sented by the wave function '4>oq---Up)o---{q')^ describes the 
configuration in which only the /i-th dressed oscillator is 
in the first excited level, all others being in their ground 
states. It is shown in Ref. [ll|, that the time evolution 
of the state \T'^) is given by 
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with 1/^1^(01 = Ij a-ll A*- This allows to inter- 
pret the coefficients f^v{t) as probability amplitudes; for 



example, fx\(t) is the probability amplitude that, if the 
dressed atom is in the first excited state at t = 0, it re- 
mains excited at time t, while f\k{t) represents the prob- 
ability amplitude that the fc-th dressed harmonic mode 
of the field be at the first excited level. 



III. TIME EVOLUTION OF A DRESSED 
TWO-ATOM STATE 

We now consider a bipartite system composed of two 
subsystems, A and B\ the subsystems consist respectively 
of dressed atoms A and i?, in the sense defined in the 
preceding section, with \ — A,B labeling the quantities 
referring to the subsystems. The whole system is con- 
tained in a perfectly refiecting sphere of radius R. In the 
following we consider each atom carrying its own dress- 
ing field (a "cloud" of field quanta), independently of 
each other. This means that we are taking the approxi- 
mation of neglecting the interaction (via the field clouds) 
between them. We consider the Hilbert space spanned 
by the dressed Fock-like product states. 



nAkik2---\ nBqiq2- 
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kiM,...) ® \^nB,qi,q2,--} ' 



I nA-ki 



(19) 



in which the dressed atom A is at the ua excited level 
and the atom B is at the tib excited level; the (doubled) 
modes of the field dressing the atoms A and B are at the 
ki,k2, ■ ■ ■ and gi, 52, • . • excited levels, respectively. Fock 
states of each individual dressed atom, A oy B, possess 
the representation and properties presented in the last 
section. 

Although it is spanned by direct products of Fock 
states of the parts, the Hilbert space of a bipartite system 
is not simply the direct product of the Hilbert spaces of 
the separated parts; it incorporates the entangled states 
as well. This is because quantum mechanics relies on the 
assumption that a linear combination of possible states 
of a given system is also an acceptable state of the sys- 
tem. Therefore, many states of a bipartite system are 
not separable, they cannot be reduced to an element of 
the direct product of the Hilbert spaces of the separated 
parts; they are entangled states which can only be con- 
ceived in a quantum mechanical framework. We shall 
now concentrate in a simple family of entangled states of 
the two dressed atom system. 

Let us consider at time t = 0, a family of superposed 
states of the bipartite system given by 
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0(A)00---; 1(B)00- 
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|U,0,0,---;0i3,0,0,---) 
+ v/r^e"^ |0^,0,0,---; 1b,0,0,---) 



(20) 
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where < ^ < 1. In this expression, |ri(^)o(s)oo - ('-')) 
and |ro'('^)i(B)oo---(0)) stand respectively for the states in 
which the dressed atom A (B) is at the first level, the 
dressed atom B (A) and all the field modes being in the 
ground state. They are 

rlS(i3)oo...(0)) = |r^oo...(o)) ® |r?oo...(0)) , (21) 
r^A?i(B)oo...(o)) = |r^oo..>|rfoo...(o)). (22) 

Note that, for ^ = 1/2 and cf) = 0,7r, states (pn| are 
similar to states of the Bell basis of a bipartite system. 

The two atoms are nondirectly interacting, they carry 
their own dressing fields (a cloud of field quanta). The 
central point, which is in the heart of the notion of en- 
tanglement, is that they share the same common wave- 
function l^f), the superposed state. In other words, we 
attribute physical reality to the superposition of the two- 
atom state I -^Q^yi^ ]^ ^^^QQ 

) , in which atom B is at the first 
excited level and the atom A in the ground state, with 
the other state Irg'^^^^'j^^^-^QQ...), in which the atom A is 
at the first excited level and the atom B in the ground 
state; afterwards, we study the time evolution of the sys- 
tem initially described by the wavefunction Eq. (PH)) . The 
field modes are all taken to be in the ground state, which 
means that we are considering the system at zero tem- 



perature. Since there is no interaction between them, the 
atoms cannot, in both classical or field-theoretical sense, 
influence one another, but as they are described by the 
same wavefunction, they are in the same superposed state 
and they can share information (not mediated by field 
forces). As largely stated in the literature, this is one of 
the more intriguing aspects of quantum mechanics; the 
correlations predicted by the theory are not compatible 
with the current idea that the state of a system, in par- 
ticular exchange of information among its subsystems, 
should be mediated by interactions among them. This 
leads still nowadays to different, yet controversial inter- 
pretations of quantum mechanics. 

In spite of the simplicity of the model, it is widely 
assumed that a pair of harmonic oscillators is a good ap- 
proximation in the case of simple atoms, for applications 
in quantum computing and for experiments with trapped 
ions. Indeed, in the realm of quantum computation i27i], 
a situation nearly equivalent to the one we investigate 
here is studied. Two noninter acting qubits, initially pre- 
pared in an entangled state, are coupled to their own in- 
dependent environments and evolve under their influence. 
This is quite similar to our approach, in which the time 
evolution of the dressed atoms is described by Eq. p7)) . 

At time t, the state of the system is described by 
the density matrix g{t) = \'^(t)) {'^{t) \ , which, using 
Eq. (PU)) . is given by 



m = e(|rfoo...W)(rfoo...W|)®(|ro^oo...)(ro^oo...|) 

+ ( 1 - g) (|r ^oo...) (r^oo...|) ® (|rfoo...(t)) {rfo,...{t)\) 

+V^(i^e^^ (|r;foo...) (rfoo...W|) ® {Ko-M (ro1,o...|) 

(|rfoo... W) (ro1)o...|) ® (Ko-) (rfoo... W|) ; (23) 



in Eq. (|23|) the states |rQQQ...) , |r^Q...) are stationary 
and the states |r^Qo... (i)), |rfoo - (i)) evolve according to 
Eq. (Ull). 

In order to investigate how the superposed states 
evolve in time, we shall consider the reduced density ma- 
trix obtained by tracing over all the degrees of freedom 
associated with the field. The computation generalizes 
the one presented in Ref . [l^] ■ After some long but rather 
straightforward calculations, we obtain the following non- 
vanishing elements 

Po^OhW = l-afAA{t)\'-{l^O\fBB{t)\\ 
PZlli*) = i^-0\fBB{t)\\ 

Phtit) = ^\fAA{t)\\ (24) 

plfiiit) - VIO^y^rAAit)fBBit), 
PifoU*) = VKl^)e~''^fAAit)fhBit)- 



We check immediately that the trace of this reduced den- 
sity matrix is one, 

pZolit)+PZll(t)+Plfolit)+Plillit) = 1: (25) 

thereby ensuring that pify represents physical states of 
the system. Also, we see that Tr [/9^(i)] ^ 1 and, there- 
fore, the superposed state at time t is not pure. The 
degree of impurity of a quantum state can be quantified 
by the departure from the idempotency property. In the 
present case, 

D(t,i) = I-Tr[p2] 

= 2(ei/AA(t)i' + (i-e)i/BB(t)i') 

-2 [^\fAA{t)\' + {l-0\fBB{t)\'y . 

(26) 

In the remainder of this section we consider the two 
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atoms as identical and, accordingly, we adopt the sub- 
script for both of them, X = A = B = 0; we also take 



fAA{t) = fBB{t) = fm{t). (27) 



9A = 9B = g; 



In this case, the matrix elements in Eqs. (|24p simplify 
and, from Eq. (I26|) . we see that the degree of impurity 
becomes independent of the superposition parameter ^: 

i?(i,0 = 2|/oo(t)|'(l-|/oo(t)|'). (28) 

In order to pursue the study of the time evolution of the 
superposition of the two-atom states, we have to deter- 
mine the behavior of foo{t). We shall analyze it in the 
situations of a very large cavity (free space) and of a small 
one. 



A. The limit of an arbitrarily large cavity 

We start from the matrix element ij^^ in Eq. (fTTl) and 
consider an arbitrarily large radius R for the cavity. The 
two (identical) atoms behave independently from each 
other, so let us focus on just one of them, either the 
atom A or the atom B, labeled 0, so that w e put A — 
A = B = 0. Remembering that rj — y^Agc/R, we have 



lim 

7?,— >oo 



lim 



y/Ag/nny/nc/R 



w2)Wl]2 



(29) 



In this limit, Auj — ttc/R dw — dQ and the sum in 
the definition of foo{t), Eq. (|18p . becomes an integral, so 
that 

/oo(t) - / dn-——-^——. (30) 
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Next, we define a parameter k — ^ Hp- ^ g^ and con- 
sider whether > q or < , for which ^ g and 
k2 -c correspond respectively to weak (g <g; w^) and 
strong ((7 3> (Da) coupling of the atoms with the envi- 
ronmment. For definiteness we consider in the following 

> 0, which includes the weak-coupling regime. We get 
in this case [lOl 



/oo(i) = e 



cos Kt — — sin Kt 



+ iG{t-Cj,g), (31) 



where the function G{t;Lj,g) is given by 

a;2 sin xt 

- I dx 

TT 



G{t-Cj,g) 



(a;2 — a;2) -|- /^g'^x'^ 

For large times, the quantity |/oo(i)|^ is given by [I^ 

2 64g2 



(32) 



i/oo(i)r 



cosut — sinwt 



(33) 



As i — !> oo, we see that the expression for |/oo(OI 8° to 
zero. 

B. Small cavity 

For a finite (small) cavity, the spectrum of eigenfre- 
quencies is discrete. Aw is large, and so the approxima- 
tion made in the case of a large cavity does not apply. 



For a sufficiently small cavity, the frequencies VLr can be 
determined as follows: in Fig. [U Eq. is plotted for 
representative values of the radius of the cavity and of 
the coupling constant. We see that apart from the small- 
est of the eigenfrequencies, all other ones are very close 
to asymptotes of the cotangent curve, which correspond 
to the field frequencies. Thus let us label the eigenfre- 
quencies as ilo, {fife}, k = 1,2,..., where fio stands to 
the smallest one. 




FIG. 1: Solutions of Eq. ((35]), with y = cot(a;) and x = 1x0.5/ g, 
for cavities satisfying the condition 5 <C 1. The asymptotes of 
the cotangent curve correspond to the frequencies of the field 
modes ujk- 



Then, defining the dimensionless parameter 
S = J— 

Aw TTC ' 

we rewrite Eq. (jl2p in the form 



cot I = — + 1 - 



g J T^g ttSU, \ g2 



(34) 



(35) 



Taking (5^1, which corresponds to i? ^ ttc/^ (a small 
cavity), we find that, for A; = 1, 2, . . ., the solutions are 



q { 26 

\ TTK 



(36) 



If we further assume that ^oirS/g <C 1, a condition com- 
patible with 5^1, then fio is found to be 



r^o ~ w 1 



(37) 



To determine /oo(i), we have to calculate the square of 
the matrix elements (^q) and . They are given, to 
first order in 5, by 



+ i^r^^'-ii^or- (38) 



We thus obtain, for sufficiently small cavities {5 ^ 1), 
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i/oowr 



1^6 



fe=l 



COS 



- ( 1 '^^ 



16^2 



7r2 ^ 



(39) 



To order 5"^, a lower bound for |/oo(i)| is obtained by 
taking the value — 1 for both cosines in the above formula, 
using the tabulated value of the Riemann zeta function 
C(2) = 7rV6: 



1 - 



47r(5 47r2,52 



(40) 



We see, comparing Eqs. ((39)) and (|33)) that the quan- 
tity |/oo(i)|^, which dictates the behavior of the density 
matrix elements and of the measure of purity in Eq. (1281) , 
has very different behaviors for free space or for a small 
cavity. This implies that in the situation of a small cav- 
ity, in contrast to the free space case, all matrix elements 
in Eqs. (p4|) are different from zero at all times. 

In Fig. ^ the degree of impurity from Ec^. (P5| is 
plotted as a function of time in the cases of an arbitrarily 
large cavity {R oo) and of a small cavity. We take 
5 = 0.1, with = 1.0 and g — 0.5 fixed (in arbitrary 
units). 




FIG. 2: Behavior of the degree of impurity D as function of 
time, equation (|28|) . for a small cavity (dashed line) and a 
very large cavity (solid line); we take the parameters g = 0.5, 
5 = 0.1 and uj = 1.0 (in arbitrary units). 



We see from the figure that for a very large cavity 
(free space) the two-atom system dissipates; with the 
passing of time, both atoms go to their ground states, 
only the element Po^o^ {t) ~ 1 survives in this limit. 
On the other hand, for a small cavity, the system never 
completely decays. 

C. Time evolution of the von Neumann entropy 

We now turn our attention to the von Neumann en- 
tropy associated with the reduced density matrix with 
respect to one of the subsystems; it is obtained by taking 
the trace over the states of the complementary subsys- 
tem in the full density matrix. For pure states of bipartite 
systems, it measures the degree of entanglement. 

The reduced density matrix for the t = Q superposition 
of states in Eq. (PH)) . p^, is obtained by tracing over the 
dressed B atom. For t ^ 0, we have 



PA{t) = Tvrs{\^{t)){-^m) 



lOO--- 



-^(l-0|^^foo...)(^^?oo...|• (4i) 



As time goes on, we have the time-dependent von Neu- 
mann entropy given by 



— — a{t) In a{t)^ 



(42) 



where here a{t) are the time-dependent eigenvalues of 
the reduced density matrix. These should be solutions 
of the characteristic equation, which in the case of (|41l) , 
reads 



det 



/1-C-a 

S.fAAfX2 









S.fA2f*AA 



UaiIX^ ei/A2|'-a 



(43) 
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We thus find that the only nonzero eigenvalues of pj^ are 

ai = l-?, a2=C^|/AMWl'=?, (44) 
/J 

which are time independent. This then implies that the 
von Neumann entropy takes the expression 

E{t, = - [(1 - i) ln(l - + e In (0] , (45) 

which coincides with the von Neumann entropy associ- 
ated with the initial state |^'(0)) given by Eq. (f20| : that 
is, all the time dependence of the von Neumann entropy 
for this two-atom system, coming from the quantities 
f\u{t), is completely cancelled in the computation of the 
entropy, in all situations, with the maximum entangle- 
ment occuring at ^ = 1/2. In other words, although the 
superposition of states evolves in time, in very different 
ways in the limits of a very large cavity and of a small 
one, the entangled nature of these two-atom states re- 
mains unchanged for all times, independently of the size 
of the cavity. 

IV. CONCLUDING REMARKS 

In this paper we have considered a system composed 
of two atoms in a spherical cavity, each of them in in- 
dependent interaction with an environment field. The 
model employed is of a bipartite system, in which each 
subsystem consists of one of the atoms dressed by its own 
proper field. We make the assumption that initially we 
have a superposition of two states: one in which one of 
the dressed atoms is in its first excited level and the other 
atom and the field modes are all in the ground state; this 
state is superposed with another one in which the atoms 
have their roles reversed. 

The time evolution of the superposition of these atomic 
states leads to a time-dependent (reduced) density ma- 
trix. Expressions for its elements are provided in both 
the cases of an infinitely large cavity (that is, free space) 
and of a small one, when the two atoms are considered 
as identical. Very different behaviors are obtained for 
this time evolution. In the large-cavity case, the system 
shows dissipation, and, with the passing of time, both 
atoms go to their ground states. For a small cavity, an 
oscillating behavior is present, so that the atoms never 
fully decay. 

In spite of these rather contrasting behaviors and of 
the nontrivial time dependence of the density matrix, we 
obtain a von Neumann entropy which is independent of 
time and of the cavity size. We find that the initial entan- 
glement of the two atoms remains unchanged as the sys- 
tem evolves, for a cavity of any size, in the approximation 
of noninteracting dressed atoms. This could be related 
to the fact that for multipartite systems the superposi- 
tion principle leads naturally to entangled states; in this 
case noninteracting subsystems can thus share entangled 
states that hold quantum correlations. Such quantum 



entanglement carries nonlocal features which can be an- 
alyzed by comparison with classical correlations [1^ lU . 
If an interaction between the dressed atoms, mediated by 
their dressing clouds, is introduced, we expect that the 
von Neumann entropy associated to the dressed atoms 
can depend on time and on the size of the cavity. How- 
ever, to establish the formalism of dressed coordinates 
and dressed states for a system of two interacting dressed 
atoms is a very hard task, which is perhaps not possible 
on purely analytical grounds. We can think of introduc- 
ing this interaction as a kind of "perturbation" around 
the individually dressed atomic states. This will be the 
subject of future work. 

We would like to emphasize that we here consider en- 
tanglement as a pure quantum effect, a characteristic of 
quantum mechanics, which is also nonlocal, in the sense 
that distant and non-interacting systems may be entan- 
gled. This is due to the existence of superposed states, 
not to the interaction between the (in our case, dressed) 
atoms. Indeed such properties of entanglement of non 
interacting systems have been used to conceive quantum 
communication devices [25| . 

Noninteracting systems have been, and currently are, 
the subject of intense investigation in the realm of tele- 
portation and quantum information theory. In (26j . en- 
tanglement in a mesoscopic structure consisting of non- 
interacting parts is investigated. These authors study 
the time-dependent electron-electron and electron-hole 
correlations in a mesoscopic device and analyze the ap- 
pearance of entanglement by means of a Bell inequal- 
ity test and of Bell inequality tests based on coincidence 
probabilities. As we have already mentioned before, in 
the framework of the theory of quantum computating, 
a situation conceptually near to the one we investigate 
here is studied [27j: two noninteracting qubits, initially 
prepared in an entangled state, are coupled to their own 
independent environments and evolve under their influ- 
ence. These authors find conditions for nonvanishing en- 
tanglement at arbitrary time, for both zero and nonzero 
temperatures. Also, in Ref. [22|, a study of the entangle- 
ment evolution of two remote atoms interacting indepen- 
dently with a cavity field is presented. In 28], quantum 
entanglement is approached for an ensemble of noninter- 
acting electrons. This author uses this as a standpoint to 
study the interacting gas and claims that in this context 
the quantum Hall effect can be thought of as a basis for 
quantum computation. 

The study of entangled states of noninteracting 
systems is interesting in itself. As clearly exposed 
in [29|], entanglement can exist as a purely quantum 
phenomenon among noninteracting particles, which 
are however described by the same wavefunction. 
Entanglement means that individual particles are not 
independent of each other, even if they do not interact, 
and their quantum properties are inextricably "tied 
up" , this being the origin of the Schrodinger's original 
denomination, verschrdnkte Zustdnde, for these states. 
In this context, the influence of an atom on the other 
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one is not due to an interaction between them, but is 
due to the attribution of physical meaning to superposed 
states, a concept with no correspondence in classical 
physics. 
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